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! Abstract 

Open charm production during the equilibration of a gluon dominated 



parton plasma is calculated, with both the time-dependent temperature and 



^ ' parton densities given by a set of rate equations. Including pre-thermal pro- 



duction, the total enhancement of open charm production over the initial 
gluon fusion depends sensitively on the initial parton density and the effec- 
tive temperature. The dependence of the pre-thermal charm production on 
the space- momentum correlation in the initial parton phase-space distribution 
is also discussed. 
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I. INTRODUCTION 



At extremely high energies, nucleus-nucleus collisions may be described by parton in- 
teractions in the framework of perturbative QCD (pQCD)-inspired models In this 
framework, hard or semihard scatterings among partons dominate the reaction dynamics. 
They can liberate partons from the individual confining nucleons, thus producing large 
amount of transverse energy in the central region |^,|^ , and drive the initially produced par- 
ton system toward equilibrium In principle, the same kind of hard processes, such 
as open charm production P|-pT|, direct photon and dilepton production [T^JT^, can also 
be used as direct probes of the early parton dynamics and the evolution of the quark-gluon 
plasma. Unlike strange quarks, charm quarks cannot be easily produced during the mixed 
and hadronic phases of the dense matter since the charm mass is much larger than the cor- 
responding temperature scale. The only period when charm quarks can be easily produced 
is during the early stage of the parton evolution when the effective temperature is still high. 
At this stage, the parton gas is still not fully equilibrated yet so that the temperature is 
only an effective parameter describing the average momentum scale. By measuring this 
pre-equilibrium charm production, one can thus probe the initial parton density in phase 
space and shed light on the equilibration time [Q. 

Roughly speaking, ultrarelativistic heavy-ion collisions in a partonic picture can be di- 
vided into three stages: (1) During the early stage, hard or semihard parton scatterings, 
which happen on a time scale of about 0.2 fm/c, produce a hot and dilute parton gas. 
This parton gas is dominated by gluons and its quark content is far below the chemical 
equilibrium value. Multiple hard scatterings suffered by a single parton during this short 
period of time when the beam partons pass through each other are suppressed due to the 
interference embedded in the Glauber formula for multiple scatterings [|1^]. This leads to 
the observed disapparence of the Cronin effect at high energy and at large transverse mo- 
mentum ll^]. Interference and parton fusion also lead to the depletion of small-x partons in 
the effective parton distributions inside a nucleus [l^JT^ . This nuclear shadowing of parton 
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distributions reduces the initial parton production |T^. (2) After two beams of partons pass 
through each other, the produced parton gas in the central rapidity region starts its evolu- 
tion toward (kinetic) thermalization mainly through elastic scatterings and expansion. The 
kinematic separation of partons through free-streaming gives an estimate of the time scale 
Tiso ~ 0.5 — 0.7 fm/c 0,0, when local isotropy in momentum distributions is reached. (3) 
Further evolution of the parton gas toward a fully (chemically) equilibrated parton plasma 
is dictated by the parton proliferation through induced radiation and gluon fusion. Due to 
the consumption of energy by the additional parton production, the effective temperature of 
the parton plasma cools down considerably faster than the ideal Bjorken's scaling solution. 
Therefore, the life time of the plasma is reduced before the temperature drops below the 
QCD phase transition temperature 0. 

Similarly, charm production can also be divided into three different contributions in the 
history of the evolution of the parton system: (1) initial production during the overlapping 
period ; (2) pre-thermal production from secondary parton scatterings during the thermal- 
ization, r < Tiso; (3) and thermal production during the parton equilibration, r > rjso, in 
the expanding system. In this paper, we will first review the equilibration of the initially 
produced parton gas in Sec. II, incorporating the result of an improved perturbative QCD 



analysis of Landau- Pomeranchuk-Migdal (LPM) effect pO| , ^ . Then we will discuss the three 
stages of open charm production in a reversed order, starting with the charm production 
during the final stage of parton equilibration in Sec. III. For pre-thermal charm produc- 
tion, we will consider the space-momentum correlation in the initial parton phase-space 
distributions, which will suppress open charm production during this period as compared 
to previous estimates ||^. In Sec. IV, we will compare the results to the charm production 
during the initial hard or semi-hard scatterings and also to the results in Geiger's calcula- 



tion |]T0[ which is about 40-50 times higher than our estimates here. We will also discuss 
the change in charm production due to the uncertainties in the initial parton density and 
effective temperature. Finally we give our conclusions in Sec. V. 
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II. PARTON PRODUCTION AND EQUILIBRATION 



At collider energies (i/i > 100 GeV), hard or semihard parton scatterings are believed 
to be the dominant mechanism for transverse energy production in the central region 0,0. 
These hard processes happen on a short time scale and they generally break color coherence 
inside individual nucleons |2^. After the fast parton pass through each other and leave the 



central region, a partonic gas will be left behind which is not immediately in thermal and 
chemical equilibrium. The partons inside such a system will then undergo further interac- 
tions and free-streaming. Neglecting parton scatterings in this period of time, the kinematic 
separation of partons with different rapidities in a cell establishes local momentum isotropy 
at the time of the order of Tiso = 0.7 fm/c |l7|Jl9(|. If we assume this is the actual kinetic 
equilibration (or thermalization) time for the partonic system. The subsequent chemical 
equilibration can then be described by a set of rate equations. In this section we will review 
parton equilibration following Ref. with improved estimate of the gluon equilibration 
rate. 

A. Initial conditions: a hot and dilute gluonic gas 

Currently there are many models for incorporating hard and semi-hard processes in 
hadronic and nuclear collisions [|T]-^. We will use the results of the HIJING Monte Carlo 
model [H] to estimate the initial parton production. In this model, multiple hard or semi- 
hard parton scatterings with initial and final state radiation are combined together with 



Lund string phenomenology p3| for the accompanying soft nonperturbative interactions. 



Let us first estimate the initial conditions at time, Tiso, from the HIJING results. Since we 
are here primarily interested in the chemical equilibration of the parton gas which has already 
reached local isotropy in momentum space, we shall assume that the parton distributions 
can be approximated by thermal phase space distributions with non-equilibrium fugacities 
A.: 
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/(fc;T,A,) = A,(e"-'=/^±A,)"\ (1) 

where is the four-velocity of the local comoving reference frame. When the parton 
fugacities Aj are much less than unity as may happen during the early evolution of the 
parton system, we can neglect the quantum corrections in Eq. (|l]) and write the momentum 
distributions in the factorized form, 

/(fc;T,A.) = A.(e-'=/^±l)"\ (2) 

Using this form of distributions, one has the parton and energy densities, 

n = {Xga, + Xgbi)T\ e = {X^a^ + Xqh2)T\ (3) 

where Oi = 16C(3)/7r^ ^ 1.95, 02 = Stt^/IS ~ 5.26, for a Bose distribution, hi = 
9C(3)A^//vr2 ^ 2.20 and 62 = 7n'^Nf/20 ^ 6.9 for a Dirac distribution. For a baryon 
symmetric system, Xq = Xg. Since boost invariance has been demonstrated to be a good 



approximation for the initially produced partons [T^ , we can then estimate the initial parton 



fugacities, A° ^ and temperature Tq from 

no = —^ 1- , €o = no-{kT}, (4) 

where (kj-) is the average transverse momentum. The quark fugacity is taken as A^ = 0.16A°, 
corresponding to a ratio 0.14 of the initial quark(antiquark) number to the total number of 
partons. Table | shows these relevant quantities at the moment Xiso, for Au + Au collisions at 
Brookhaven National Laboratory Relativistic Heavy Ion Collider (RHIC) and CERN Large 
Hadron Collider (LHC) energies. One can observe that the initial parton gas is rather hot 
reflecting the large average transverse momentum. However, the parton gas is very dilute 
as compared to the ideal gas at the same temperature. The gas is also dominated by gluons 
with its quark content far below the chemical equilibrium value. We should emphasize 
that the initial conditions listed here result from HIJING calculation of parton production 
through semihard scatterings. Soft partons, e.g., due to parton production from the color 



field [24], are not included. 



B. Master rate equations 



In general, chemical reactions among partons can be quite complicated because of the 
possibility of initial and final-state gluon radiations. Interference effects due to multiple 
scatterings inside a dense medium, i.e., LPM suppression of soft gluon radiation has to be 
taken into account. One lesson one has learned from LPM effect pO| , pl| is that the radiation 
between two successive scatterings is the sum, on the amplitude level, of both the initial 
state radiation from the first scattering and the final state radiation from the second one. 
Since the off-shell parton is space-like in the first amplitude and time-like in the second, the 
picture of a time-like parton propagating inside a medium in the parton cascade simulations 
0] shall break down. Instead, we shall here consider both initial and final state radiations 
together associated with a single scattering (To the leading order, a single additional gluon is 
radiated, such as gg —>■ ggg) , in which we can include LPM effect by a radiation suppression 



factor. The analysis of QCD LPM effect in Ref. ||20|j21| has been done for a fast parton 
traveling inside a parton plasma. We will use the results for radiations off thermal partons 
who average energy is about T, since we expect the same physics to happen. 

In order to permit the approach to chemical equilibrium, the reverse process, i.e., gluon 
absorption, has to be included as well, which is easily achieved making use of detailed 
balance. We consider only the dominant process gg —>■ ggg. Radiative processes involving 
quarks have substantially smaller cross sections in pQCD, and quarks are considerably less 
abundant than gluons in the initial phase of the chemical evolution of the parton gas. 
Here we are interested in understanding the basic mechanisms underlying the formation of a 
chemically equilibrated quark-gluon plasma, and the essential time-scales. We hence restrict 
our considerations to the dominant reaction mechanism for the equilibration of each parton 



flavor. These are just four processes [25 



99 ^ 999, 99 ^ QQ- (5) 

Other scattering processes ensure the maintenance of thermal equilibrium {gg ^ gg, gq ^ 
gq, etc.) or yield corrections to the dominant reaction rates {gq ^ qgg, etc.). 
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Restricting to the reactions in Eq. and assuming that elastic parton scatterings are 
sufficiently rapid to maintain local thermal equilibrium, the evolution of the parton densities 
is governed by the master equations M: 



Mn,un ^ (l-^^y (l - Jj) . (6) 

iUn,nn = \a,nl (l " ||) . (7) 

where hi = Ui/Xi denote the densities with unit fugacities, Aj = 1, cxs and a2 are thermally 
averaged, velocity weighted cross sections, 

(^3 = {(^{99 999)v), (^2 = W{99 (8) 

We have also assumed detailed balance and a baryon symmetric matter, Uq = Uq. If we 
neglect effects of viscosity due to elastic scattering ||2^ , |26[| , we then have the hydrodynamic 
equation 



d^{eu^) + P d^u^ = (9) 

which determines the evolution of the energy density. 

For a time scale r -C Ra-, we can neglect the transverse expansion and consider the 
expansion of the parton plasma purely longitudinal, which leads to the Bjorken's scaling 



solution of the hydrodynamic equation: 



^ + ^^ = 0. (10) 

We further assume the ultrarelativistic equation of motion, e = 3P with rii and e given 
by Eq. (^). We can then solve the hydrodynamic equation, 

[Xg + -Xgf^T^T = const. , (11) 

0,2 

and rewrite the rate equation in terms of time evolution of the parameters T(r), Ag(r) and 
\(r), 
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^ + 3| + i^fi3(l-A,)-2«.(l-|) (12) 

where the density weighted reaction rates Rs and R2 are defined as 

^3 = |o-3^s, R2 = \(y2ng. (14) 

Notice that for a fully equilibrated system (A^ = = 1), Eq. ( ]TT| ) corresponds to the 
Bjorken solution, T{t)/Tq = (tq/tY^^. 



C. Parton equilibration rates 

To take into account of the LPM effect in the calculation of the reaction rate R3 for 
99 ~^ 999i we simply impose the LPM suppression of the gluon radiation whose effective 
formation time tqcd is niuch longer than the mean-free-path A/ of multiple scatterings to 
each gg ggg process. In the mean time, the LPM effect also regularizes the infrared di- 
vergency associated with QCD radiation. However, still contains infrared singularities in 
the gluon propagators. For an equilibrium system one can in principle apply the resumma- 
tion technique developed by Braaten and Pisarski pS| to regularize the electric part of the 



propagators, though the magnetic sector still has to be determined by an unknown magnetic 
screening mass which can only be calculated nonperturbatively up to now. Since we are 
dealing with a nonequilibrium system, Braaten and Pisarski's resummation may not be well 
defined. As an approximation, we will use the Debye screening mass |2^ , 

A = % kfik)dk = ATiasT^Xg, (15) 

to regularize all singularities in both the scattering cross sections and the radiation ampli- 
tude. 

To further simplify the calculation we approximate the LPM suppression factor in 
Ref. [pn| , PT| by a 6'-function, 9{\f — tqcd)) where 
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Ca 2 coshy 

= 2C^^^' ^^^^ 

is the effective formation time of the gluon radiation in QCD which depends on the second 
Casimir of the beam parton representation in SU{3), e.g., C2 = Ca = 3 for a gluon. In the 
previous calculation of the interaction rate [0, this color factor was not taken into account. 
The modified differential cross section for gg —>■ ggg is then, 

where the second step-function accounts for energy conservation, and s = 18T^ is the average 
squared center-of-mass energy of two gluons in the thermal gas. The regularized gluon 



density distribution induced by a single scattering is |30 

dUg CaUs ql 



Similarly, the regularized small angle gg scattering cross is. 



(18) 



dql Hql+fiiy 
The mean-free-path for elastic scatterings is then. 



(19) 



A7^ = ug / dql-f- = -aia.T— — — -, (20) 

■' Jo dq^ 8 1 + 8nasXg/9 

which depends very weekly on the gluon fugacity as compared to the independent one 
used in a previous study Using, 

. , 1 27r 



(21) 



we can complete part of the integrations and have, 

32 

R3/T = —asXgil + 87ias\g/9Yli\g), (22) 



where ^(Ag) is a function of A 



9' 



9 



, /v^^/ , /-^/^Md , z \ cosh ^iJx) 

I{\g) = dx dz- > ^ 



1 + zy ^ xJ[x + (1 + z)xd]'^ — Ax z xd 



+ ^ , } , (23) 

s^f J[l + x{l + z)yDf -Ax zyo^ 



where xd = fJqj'^f ^^^d i/d = /^d/s. We can evaluate the integration numerically and find 
out the dependence of R3/T on the gluon fugacity Xg. In Fig. H R3/T is plotted versus 
for a coupling constant = 0.3. The gluon production rate increases with and then 
saturates when the system is in equilibrium. Note that in principle one should multiply the 
phase-space integral by 1/3! to take into account of the symmetrization of identical particles 
in gg ggg as in Ref. |T^. However, for the dominant soft radiation we consider here, the 
radiated soft gluon does not overlap with the two incident gluons in the phase-space. Thus 
we only multiply the cross section by 1/2! to obtain the equilibration rate. 

The calculation of the quark equilibration rate R2 for gg qq is more straightforward. 
Estimate in Ref. gives, 

R2 = ^(T2ng ^ 0.2ANf a%T\n{5.5/\g). (24) 

The dashed line in Fig. |l| shows the normalized rate R2/T for Nf = 2.5, taking into account 
the reduced phase space of strange quarks at moderate temperatures, as a function of the 
gluon fugacity. 



D. Evolution of the parton plasma 

With the parton equilibration rates which in turn depend on the parton fugacity, we can 
solve the master equations self-consistently and obtain the time evolution of the temperature 
and the fugacities. Shown in Figs. ^ and ^ are the time dependence of T, A^, and A^ for 
initial conditions listed in Table | at RHIC and LHC energies. We find that the parton gas 
cools considerably faster than predicted by Bjorken's scaling solution {T^t = const.) shown 
as dotted lines, because the production of additional partons approaching the chemical 
equilibrium state consumes an appreciable amount of energy. The accelerated cooling, in 

10 



turn, slows down the chemical equilibration process, which is more apparent at RHIC than 
at LHC energies. Therefore, the parton system can hardly reach its equilibrium state before 
the effective temperature drops below Tc ~ 200 MeV in a short period of time of 1-2 fm/c at 
RHIC energy. At LHC energy, however, the parton gas becomes very close to its equilibrium 
and the plasma may exist in a deconfined phase for as long as 4-5 fm/c. Another important 
observation is that quarks never approach to chemical equilibrium at both energies. This 
is partially due to the small initial quark fugacity and partially due to the small quark 
equilibration rate. 

We note that the initial conditions used here result from the HIJING model calculation 
in which only initial direct parton scatterings are taken into account. Due to the fact 
that HIJING is a pQCD motivated phenomenological model, there are some uncertainties 
related to the initial parton production, as listed in Ref . . We can estimate the effect of the 
uncertainties in the initial conditions on the parton gas evolution by multiplying the initial 
energy and parton number densities at RHIC energy by a factor of 4. This will result in 
the initial fugacities, A° = 0.2 and A° = 0.024. With these high initial densities, the parton 
plasma can evolve into a nearly equilibrated gluon gas as shown in Fig. ^. The deconfined 
phase will also last longer for about 4 fm/c. Though, the system is still dominated by gluons 
with few quarks and antiquarks as compared to a fully chemical equilibrated system. If the 
uncertainties in the initial conditions are caused by the soft parton production from the color 
mean fields, the initial effective temperature will decrease. Therefore, we can alternatively 
increase the initial parton density by a factor of 4 and decrease Tq to 0.4 GeV at the same 
time. This leads to higher initial fugacities, A° = 0.52 and Ajj = 0.083. As shown in Fig. ^by 
the curves with stars, this system evolves faster toward equilibrium, however, with shorter 
life-time in the deconfined phase due to the reduced initial temperature. 

We thus can conclude that perturbative parton production and scatterings are very likely 
to produce a quark-gluon plasma ( or more specifically a gluon plasma) in ultrarelativistic 
heavy ion collisions at LHC energy. However, the fate of the quark-gluon plasma at RHIC 
energy has to be determined by a more careful examination of the uncertainties in the 
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initial conditions. These uncertainties will surely affect the open charm production during 
the equilibration as we shall discuss. 



III. THERMAL CHARM PRODUCTION DURING EQUILIBRATION 



With the given evolution of the parton gas, we can now calculate open charm production 
during the parton equilibration. Similar to light quarks, charm quarks are produced through 
gluon fusion gg — > cc and quark antiquark annihilation qq cc during the evolution of the 
parton plasma. However, since the number of charm quarks is very small as compared to 
gluons and light quarks, we can neglect the back reactions, cc gg, qq and their effect on 
the parton evolution. Given the phase-space density of the equilibrating partons, fi{k), the 
differential production rate is 0, 



E- 



(Pp 16(27r)8 
1 



(fki <fk2 d?P2 

LOl UJ2 E2 



6^^\h + k2-p- 



P2) 



i:9ygih)fg{k2)\Mgg^Cc\^ + g^Jg{kl)f^{k2)\Mgg-,cc\ 



(25) 



where (^^=16, gq = 6Nf, are the degeneracy factors for gluons and quarks (antiquarks) 
respectively, {Aigg^ccl"^, lAigq^ccl"^ are the averaged matrix elements for gg cc and qq cc 
processes, respectively, 

\Mgn^rr\'' VI,.,. 8 f M'' - U - t 



gg^,^,^ ^ 12 



u 



8 (M^ -u M2 
3 \M^-i^ M^-u, 



16M2 



-(2M2 

s 



i-u) 



6 

+ - 



M\i-uf 



2 Af^is-AM^) 



s (M2 - t) (M2 - u) 3 (M2 - t) (M2 - u) 



\M 



qq^cc\ 



64 
9P 



{M^ -if + {M^ -iif + 2M^s 



(26) 
(27) 



Due to the small charm density, we can neglect the Pauli blocking of the final charm quarks. 
For large charm quark mass, M, we can approximate the phase-space density fi{k) by a 
Boltzmann distribution. We further assume that the distributions are boost invariant, i.e., 



(28) 
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where rj — 0.5 ln{t + z)/{t — z) is the spatial rapidity of a space-time cell at {t,z). Neglecting 
the transverse expansion, the above assumption implies that the space-time cell at (t, z) 
have a flow velocity, u = (coshr;, sinhr;). We can now complete the integral over 77 in 
/ d^x = 7ri?^ / drjdT and obtain, 

^ I I P±2dp±2d(f)2dy2d(f)kidyki'^^^2Ko{Q±/T) 

Lb(^7rj° Jn^ J s 

■1 „ „i ^^^^ 



where Kq is the modified Bassel function and Tc is the time when the temperature, T, drops 
below 200 MeV. The kinematic variables are chosen such that. 



P2 = (Mj_2 cosh 1/2, p±2 COS (/)2, p±2 siu 02, Mj_2 siuh 112), Mj_2 = y'M^ p^2; 

ki = /i;±i(cosh Uki, cos sin sinh yki), i = 1, 2 . (30) 

The center-of-mass momentum, Q = {Q± coshyg, q^, Q± sinhyg), is defined as Q = p+P2 — 
ki + k2, and 

q2 ^ s = 2[M2 + MrM^2 cosh(7/ - ^2) - P±P±2 cos ^2], 

q1^P±+ PI2 + 2p±P±2 cos (f)2, 

Ql = Q^ + ql = M^ + Ml^ + 2MtM^2 cosh(y - 7/2)- (31) 
Using these variables and the energy-momentum conservation, we have, 

gV2 



k 



±1 



Mj_ cosh(|/ - yki) + Mj_2 cosh(y2 - Vki) - q± cos (pig ' 
cos 0ig = [p_L cos + P±2 cos(02 - 4>ki)]/q±- (32) 

In the integral over r, we shall use the time evolution of the temperature, T(r), and fugac- 
ities, Ai(T), as given in the previous section. 



IV. PRE-THERMAL CHARM PRODUCTION 



Before the parton distributions reach local isotropy in momentum space so that the rate 
equations can be applied to describe the equilibration of the parton system, scatterings 
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among free-streaming partons can also lead to charm production. Since the system during 
this period consists dominantly of gluons, we shall only consider gluon fusions. To model 
the phase-space distribution, we take into account the distribution of the initial production 
points which spread over a region of width, 



in z coordinate. Following Ref. |TT[, we assume free-streaming until t-^so and neglect the 
expansion in the transverse direction. The correlated phase-space distribution function is 
given by 

2 g-(2:-ttanhy)2/2A2 

f{k,x) = ^g{k±,y) j==- e{RA - r)e{Ti,o cosh y - t), (34) 

where g{k±,y) is the parametrization of the parton spectrum given by HIJING simulations, 

The phase-space distribution is normalized such that lim^^oo S'g / '^'^3;/(/c, x)/(27r)^ = 
d?Ng/d^k. The function h{k±) and the rapidity width Y are given in Table |T| for cen- 
tral Au + Au collisions at RHIC and LHC energies which also gives the initial conditions as 
listed in Table |. 

Substituting the phase-space distribution, into Eq. (p5|) , and integrate over space and 
time, we obtain the charm production distribution in the pre-thermal period. 



diV if Zk 

dyd^p_[_ 16(27r)"7ritji J s 

2 \/27rA+„+ Jo 



Hot - \/A|i + AI2, tf = Tiso min(cosh?/fei, cosh?/fc2)- (37) 
where the kinematic variables are similarly defined as in Eq. (p9D, and in addition, 

^±2 = Q± + kli - 2/c_Li[M_L cosh(?/ - yki) + Mj_2 cosh(?/2 - 
sinh yk2 = [M± sinh y + M±2 sinh 2/2 - ^±1 sinh yki]/k±2. (38) 
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Note that the correlation between momentum and space-time in the phase-space distri- 
bution was not considered in a previous calculation 0. As we will show this correlation 
is very important and will reduce the pre-thermal charm production as compared to the 
uncorrelated distributions. Similar effect was recently discussed by Lin and Gyulassy in 



Ref . [O , where formation time effect is also included which is expected to further suppress 



pre-thermal charm production. 

V. INITIAL FUSION 

During the initial interaction period, charm quarks are produced together with minijets 
through gluon fusion and quark anti-quark annihilation. Like gluon and light quark produc- 
tion, charm production through the initial fusion is very sensitive to the parton distributions 
inside nuclei. In addition, the cross section is also very sensitive to the value of charm quark 
mass, M. If higher order corrections are taken into account, the production cross section 
depends also on the choices of the renormalization and factorization scales. Detailed studies 



on the next-leading-order calculation [pl| -|3^ shows, however, that higher order corrections 
to the total charm production cross section can be accounted for by a constant A"-factor 
of about 2. This is what we will use next. For consistency we use M = 1.5 GeV for all 
calculations. Shown in Figs. |^ and ^ as solid lines are the initial charm production given 
by HIJING calculations at RHIC and LHC energies, with MRSD— ' parton distribu- 
tions. The corresponding total integrated cross sections are, Ucc = 0.16 (5.75) mb at RHIC 
(LHC) energy, where nuclear shadowing of the gluon distribution function is also taken into 
account. In HIJING calculations, high order corrections are included via parton cascade in 
both initial and final state radiations. The resultant distributions in cc-pair momentum are 
very close to the explicit higher order calculations [p4| . 

Plotted in Figs. |^ and ^ as dot-dashed and dashed lines are the pre-thermal and thermal 
production. In the calculation, a factor of 2 is also multiplied to the lowest order matrix 
elements of charm production. Both contributions are much smaller than the initial charm 
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production at both energies. The pre-thermal contributions shown are also much smaller 
than what was found in Ref. 0. This is because momentum and space-time correlation 
was not taken into account in Ref. ^ which suppresses the pre-thermal charm production. 
Similar results are also found in a study by Lin and Gyulassy As we have already 

discussed, the initial conditions in Tables |I| and ^ given by HIJING calculations have many 
uncertainties. If one increases the initial parton number density at RHIC energy by a 
factor of 4 with the same initial temperature, charm production from both pre-thermal and 
thermal sources will increase about a factor of 12 as shown in Fig. 0, leading to a total 
secondary contribution comparable to the initial charm production. In the extreme limit, 
a fully equilibrated parton plasma (A^ = Ag = 1) at the same initial temperature would 
give an enhancement of charm production about 4 times higher than the initial production, 
shown as dotted lines in Figs. ^ and |^. In this case, the enhancement not only comes from 
higher parton densities, but also from the much longer life time of the parton plasma (cf. 
Figs. and 0). Much higher enhancements predicted in Ref. |T0| are due to the overestimate 



of the intrinsic charm production as pointed by Gyulassy and Lin in a recent paper [IT 
Though the intrinsic charm production is important in the forward direction at large x j 



it is strongly suppressed in the mid-rapidity region due to the interference among pQCD 



amplitudes to the same order |37 



To test the sensitivity of open charm production to uncertainties in initial fugacities 
and temperature separately, we consider an alternative scenario as we have discussed in the 
parton evolution. We assume the initial parton densities to be 4 times higher than given 
in Table |I| at RHIC energy but with lower initial temperature, Tq = 0.4 GeV. Accordingly, 
the initial phase-space distribution is also modified to: h{k±) = 9649.2e~^^/°-^V(A;x + 0.3) 
from the one in Table y, which gives 4 times of the initial parton density but smaller av- 
erage transverse momentum, = 0.85 GeV. The reduced average transverse momentum 
corresponds to lower initial effective temperature. This system with higher initial fugacities 
evolves faster toward equilibrium but the life-time of the deconfined phase is shorter due to 
the reduced temperature as we have discussed. The corresponding open charm production 

16 



is shown in Fig. |^ by the hnes with stars. We observe that open charm production from 
both pre-thermal and thermal contribution is reduced due to the reduction in initial temper- 
ature and life-time of the parton plasma, even though the initial fugacities are much higher 
and the evolution toward equilibrium is faster. Thus, open charm production is much more 
sensitive to the change in the initial temperature than the parton fugacities. We also note 
from Eqs. (p9|) and (|36|) that the pre-thermal and thermal charm production depends on the 
thermalization time tiso and the life time of the parton plasma. Therefore, by measuring 
the charm enhancement, we can probe the initial parton phase-space distribution, initial 
temperature and the thermalization and equilibration time. 



VI. CONCLUSIONS 

In this paper, we have calculated open charm production in an equilibrating parton 
plasma, taking into account the evolution of the effective temperature and parton fugacities 
according to the solution of a set of rate equations. In the evaluation of the interaction 
rate R3 for induced gluon radiation, a color dependent effective formation time was used 
which reduces the gluon equilibration rate through LPM suppression of soft gluons. In the 
calculation of the pre-thermal contribution to open charm production, correlation between 
momentum and space-time was also included. This correlation reduces the pre-thermal 
charm production as compared to the uncorrelated one used in a previous estimate 0. 

We found that both the thermal contribution during the parton equilibration and pre- 
thermal contribution with the current estimate of the initial parton density from HIJING 
Monte Carlo simulation are much smaller than the initial direct charm production. However, 
the final total charm production is very sensitive to the initial condition of the parton 
evolution. If uncertainties in the initial parton production can increase the initial parton 
density, e.g., by a factor of 4, the total secondary charm production will become comparable 
or larger than the initial production, due to both the increased production rate and longer 
life time of the parton plasma. We also found that open charm production is more sensitive 
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to the initial temperature of the parton system than the initial parton fugacities. Therefore, 
open charm production is a good probe of the initial parton distribution in phase-space and 
the thermalizaion and equilibration time of the parton plasma. 
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TABLES 



EIIIC LIIC 

Tiso (fm/c) 0.7 0.5 

£o (GeV/fm3) 3.2 40 

no (fm-=^) 2.15 18 

{k±) (GeV) 1.17 1.76 

To (GeV) 0.55 0.82 

A° 0.05 0.124 

A° 0.008 0.02 

TABLE I. Values of the relevant parameters characterizing the parton plasma at the moment 
Tiso, when local isotropy of the momentum distribution is first reached. 



^/s (TeV) Y h{k±_) (GeV- 2) 

0.2 2.5 1754.4e-'=^/°-7(fe± + 0.3) 

5.5 5.0 2.66 X 107(A;± + 2.9)^-^ 

TABLE IL Parametrizations of the momentum spectra of the initially produced partons in 
HIJING calculation. 
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FIGURES 

FIG. 1. The scaled gluon production rate R3/T (solid line) for gg ggg and the quark 
production rate R2/T (dashed line) for gg — > qq are shown as function of the gluon fugacity Xg for 
as = 0.3. 

FIG. 2. Time evolution of the temperature T and the fugacities and Xg of gluons and quarks 
in the parton plasma created in Au + Au collisions at the RHIC energy of ^/s = 200 AGeV. The 
initial values for T, Xg and Xq are determined from HIJING simulations and are listed in Table |. 

FIG. 3. The same as in Fig. ^, except for LHC energy, ^/s = 5.5 ATeV. 

FIG. 4. The same as in Fig. ^, except that the initial parton densities are 4 times higher than 
given in Table || with the same (ordinary lines), or reduced initial temperature, Tq = 0.4 GeV (lines 
with stars) 

FIG. 5. The p± distributions of initial (solid), prethermal (dot-dashed), and thermal (dashed) 
charm production for central Au + Au collisions at RHIC energy, ^/s = 200 AGeV with initial 
conditions given in Tables | and The dotted line is the thermal production assuming an initial 
fully equilibrated QGP at the same temperature. 

FIG. 6. The same as in Fig. |5[ except at LHC energy, ^/s = 5.5 ATeV. 

FIG. 7. The p± distribution of the initial (solid), pre-thermal (dot-dashed) and thermal 
(dashed) charm production for initial parton densities 4 times higher than HIJING estimate given 
in I but with the same (ordinary lines), or reduced initial temperature, Tq = 0.4 GeV (lines with 
stars). 
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